Nuclear matter fourth-order symmetry energy in the relativistic mean field model 
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Within the nonhnear relativistic mean field model, we derive the analytical expression of the 
nuclear matter fourth-order symmetry energy -Esym,4(p)- Based on two accurately calibrated in- 
teractions FSUGold and lU-FSU, our results show that the value of -Bsym,4(p) at normal nuclear 
matter density po is generally less than 1 MeV, confirming the empirical parabolic approximation 
to the equation of state for asymmetric nuclear matter at po- On the other hand, we find that 
the iJsym,4(p) may become nonnegligible at high densities. Furthermore, the analytical form of the 
£^sym,4(p) provides the possibility to study the higher-order effects on the isobaric incompressibility 
of asymmetric nuclear matter, i.e., KaatiS) — Ko + KBa.t,2S^ + KsB.t,45* + 0{6^) where S = ipn — pp)/p 
is the isospin asymmetry, and we find that the value of -R'sat,4 is generally small compared with that 
of the -R'sat,2. In addition, we study the effects of the i5syni,4(p) on the proton fraction Xp and 
the core-crust transition density pt and pressure Pt in neutron stars. Interestingly, we find that, 
compared with the results from the empirical parabolic approximation, including the Esyin,i{p) con- 
tribution can significantly enhance the Xp at high densities and strongly reduce the pt and Pt in 
neutron stars, demonstrating that the widely used empirical parabolic approximation may cause 
large errors in determining the Xp at high densities as well as the pt and Pt in neutron stars within 
the nonlinear relativistic mean field model, consistent with previous nonrclativistic calculations. 

PACS numbers: 21.65.Ef, 24.10.Jv, 26.60.Gj, 21.30.Fe 



I. INTRODUCTION 



One of fundamental issues in nuclear physics is the 
equation of state (EOS) of isospin asymmetric nuclear 
matter, which plays a central role in understanding not 
only the structure of radioactive nuclei, the reaction dy- 
namics induced by rare isotopes, and the liquid-gas phase 
transition in asymmetric nuclear matter, but also many 
critical issues in astrophysics [ll-Q]- For symmetric nu- 
clear matter with equal fractions of neutrons and pro- 
tons, its EOS is relatively well-determined from analyses 
of the giant monopole resonances of finite nuclei 8, 9] as 
well as collective flows and subthreshold kaon produc- 
tion [13, [U in relativistic nucleus-nucleus collisions. On 
the other hand, the EOS of asymmetric nuclear matter, 
especially the density dependence of the nuclear symme- 
try energy i?sym(p), is poorly known. During the last 
decade, significant progress has been made both experi- 
mentally and theoretically on constraining the behavior 
of the symmetry energy around and below normal nu- 
clear matter density [la - fTsj (See, e.g., Refs. [l6l - [20j for 
review of recent progress) while its super-normal density 
behavior remains elusive and largely controversial |2l| - 
[23 |. Theoretically, all many-body theory calculations to 
date have demonstrated that the nuclear symmetry en- 
ergy essentially characterizes the isospin dependent part 
of the EOS of asymmetric nuclear matter and the higher- 
order terms in isospin asymmetry are unimportant, at 
least for densities up to moderate values Q, leading to 
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the well-known empirical parabolic law. 

When the empirical parabolic law itself provides a good 
approximation to the EOS of asymmetric nuclear mat- 
ter and thus allows one to extract the symmetry energy 
from the energy difference between pure neutron matter 
and symmetric nuclear matter, it may cause large errors 
when it is applied to determine some physical quanti- 
ties under special conditions. For example, the higher- 
order terms in isospin asymmetry presented in the EOS of 
asymmetric nuclear matter at supra-normal densities can 
significantly modify the proton fraction in /3-equilibrium 
neutron-star matter and the critical density for the di- 
rect Urea process which can lead to faster cooling of neu- 
tron stars [1^, [2^. In addition, recent studies [l^l in- 
dicate that the higher-order terms in isospin asymmetry 
are very important for determining the transition den- 
sity and pressure at the inner edge separating the liquid 
core from the solid crust of neutron stars where the mat- 
ter is extremely neutron-rich. Furthermore, the higher- 
order effects on the incompressibility of asymmetric nu- 
clear matter have also been studied recently [l^. These 
studies about the higher-order effects are essentially per- 
formed within the nonrelativistic models since the ana- 
lytical expressions of the higher-order terms in isospin 
asymmetry, e.g., the nuclear matter fourth-order sym- 
metry energy £'sym,4(p), can be relatively easily obtained 
in such nonrelativistic models. It is thus interesting to 
see if the same conclusion can be obtained within the 
relativistic models. 

One of very popular relativistic models is the rel- 
ativistic mean field (RMF) model which is generally 
based on effective interaction Lagrangians involving nu- 
cleon and meson fields [2^. As a phenomenological ap- 



2 



proach, the RMF model has achieved great success dur- 
ing the last decades in describing many nuclear phenom- 
ena [13, HH. Although the full expressions have been 
usually used in realistic RMF model calculations, never- 
theless, it will be instructive to see separately the effects 
of the higher-order terms in isospin asymmetry, e.g., the 
nuclear matter fourth-order symmetry energy Esym,4:ip)j 
within the RMF model. The main motivation of the 
present work is to derive the nuclear matter fourth-order 
symmetry energy £^sym,4(p) within the nonlinear RMF 
model and then explore the higher-order i?sym,4(p) cor- 
rections to the widely used empirical parabolic law for 
the isospin asymmetric nuclear matter. Based on two 
accurately calibrated interactions, our results indicate 
that the EgymAip) may have significant influence on the 
properties of isospin asymmetric nuclear matter, the pro- 
ton fraction Xp in /3-stable neutron star matter and the 
core-crust transition density pt and pressure Pt in neu- 
tron stars, confirming the previous nonrelativistic calcu- 
lations. 

The paper is organized as follows. In Section 
we briefiy discuss the bulk characteristic parameters of 
asymmetric nuclear matter. The model used in the 
present paper and the analytical expression of the 4th- 
order symmetry energy will also given in this section. 
The results and discussions are then presented in Section 
mil Finally, a summary is given in Section IIVI 



S are usually very small. For example, the magnitude of 
the (5^ term at normal nuclear matter density po is esti- 
mated to be less than 1 MeV in microscopic many-body 
approaches |32l43JI and also in phenomenological nonrel- 
ativistic models |28| as well as relativistic models as will 
be shown in this work. Neglecting the contribution from 
higher-order terms in Eq. ([T]) leads to the well-known em- 
pirical parabolic law, i.e., E{p,6) ~ EQ{p)+Esyni{p)S^ for 
the EOS of asymmetric nuclear matter and the symme- 
try energy i?sym(p) can thus be extracted from i?sym(p) — 
E{p,S = l)-E{p,S^O). 

Around normal nuclear matter density po, the Eq(p) 
can be expanded, e.g., up to 4th-order in density, as. 



Eo{p) = Eoipo) 



fx^ + fx-V§x^ + 0(x^), (4) 



where x is a dimensionless variable characterizing the de- 
viations of the density from normal nuclear matter den- 
sity Po £md it is conventionally defined as 



X 



3po 



(5) 



The first term i?o(po) on the right-hand-side (r.h.s) of 
Eq. (ID) is the binding energy per nucleon in symmetric 
nuclear matter at normal nuclear matter density po and 
the coefficients of other terms are, 



II. THEORETICAL FORMULISM 



Characteristic parameters of asymmetric 
nuclear matter 



The EOS of isospin asymmetric nuclear matter, defined 
by its binding energy per nucleon, can be expanded to 
4nd-order in isospin asymmetry 6 as 

E{p,6) = Eoip) + Esy^{p)6^ + EsymAip)S^ + 0(6''), (1) 

where p ~ Pn + Pp is the baryon density with p„ and 
Pp denoting the neutron and proton densities, respec- 
tively; S ~ {pn — Pp)/{pp + Pn) is the isospin asymmetry; 
Eo{p) ~ E{p,5 — 0) is the binding energy per nucleon 
in symmetric nuclear matter; the nuclear matter symme- 
try energy -Bsym(p) and the 4th-order symmetry energy 
i?sym,4(p) are expressed, respectively, as 



Esymip) 
EsymAiP) 



1 d^Eip,S) 



2! dS^ 
1 d^Eip,S) 



4! dS^ 



5=0 



5=0 



(2) 
(3) 



In Eq. ((T|), the absence of odd-order terms in 5 is due to 
the exchange symmetry between protons and neutrons in 
nuclear matter when one neglects the Coulomb interac- 
tion and assumes the charge symmetry of nuclear forces. 
The higher-order (including the 4th-order) coefficients in 



A'o 
Jo 
lo 



9p, 



2d^E^{p) 



dp^ 



dp^ 



dp* 



(6) 
(7) 
(8) 



The linear x term on the r.h.s of Eq. ^ vanishes accord- 
ing to the definition of the saturation density pQ. The 
coefficient A'o is the well-known incomprcssibility coef- 
ficient of symmetric nuclear matter and it characterizes 
the curvature of Eo{po) at po- The coefficients Jo and Jo 
are the 3rd-order and 4th-order incompressibility coeffi- 
cients of symmetric nuclear matter (28j . respectively. 

Similarly, around normal nuclear matter density po, 
the symmetry energy Egy^ip) and the 4th-order symme- 
try energy EsymA{p) can be expanded, e.g., up to 4th- 
order in as, 



-E'sym (p) = 



Esym(,Po) 



Lx 



_sym 2 

2! ^ 



-^x + 



3! 



4! 



x' + 0{x'), (9) 



and 



E. 



'sym,4 



(p) 



E. 



sym,4 



(po) 



-^sym,4 2 

2! ^ 



3! ^ 



L 



sym,4 4 



4! 



x' + oixn. (10) 
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respectively, where the L, ii'sym, 



K 



sym,4; <^sym,4; -^sym,4 



sym; -'sym 



and L 



sym,4 ; 



are the slope parameter, curvature 
parameter, 3rd-order and 4th-order density coefRcients 
of the Esymip) and i?sym.4(p) at po, respectively, whose 
definitions are similar to Eq. (|6]) - Eq. (|8]). In general, 
these characteristic parameters can be written as, 



(3po) 



(11) 



P=PO 



for example, Wn — L, W\i 



W- 
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J. 



sym,4 1 



W2^ 



sym,4 5 



and so on. 



In the above Taylor's expansions, we have kept all 
terms up to 4th-order in 5 or x- The 14 characteristic pa- 
rameters, namely, £'o(po), Kt^, Jq, /q, i^sym(po), Ksym, 



J. 



sym; -^syin; -^syin,4(po): -^sym,4: -^sym,4 



; '^sym,4 and /sym,4 



are well-defined, and they characterize the EOS of an 
asymmetric nuclear matter and its density dependence at 
normal nuclear matter density pQ. Among these param- 
eters, Eo{po), Kf), Esymipa), L and Kgym have been ex- 
tensively studied in the literature and significant progress 
has been made over past few decades [7|- 

The incompressibility of asymmetric nuclear matter 
is an important quantity to characterize its EOS. Con- 
ventionally, the incompressibility coefficient is defined 
at the saturation density where the pressure P{p, S) = 
p^dE{p,6)/dp = 0, and it is called the isobaric incom- 
pressibility coefficient [s^ given by 



d^E{p,S) 



9p2 



(12) 



P=Psat 



The isobaric incompressibility coefficient A'satlf^) thus 
only depends on the isospin asymmetry 6. One can show 
that up to 4th-order in 5, the Kgati^) can be expressed 
as [H 



A'sat,4<5^ 



with 



sat, 2 



Ksym 6i 



-^sat,4 — ^sym,4 



6Z/sym,4 

JoKsymL 
^0 



JqL, 



sym,4 



^0 



9L^ 
AcT 
4L 



(13) 
(14) 

in L 



Kr 



2Kl 



(15) 



If we use the parabolic approximation for the EOS of 
symmetric nuclear matter, i.e., Eq{p) = Eo{pq) + ^KqX~^ 
0{x^), then the A'sat,2 parameter is reduced to 



B. The 4th-order symmetry energy in the 
nonlinear RMF model 



In the present work, we use the interacting Lagrangian 
density of the nonlinear RMF model supplemented with 
couplings between the isoscalar and the isovector mesons 
[1311, i.e., 

£(?/', cr,a;,,, p^) = i' [lt,{id^ ~ 9^'^^) - [M - g^a)] 



-g^p^ ■ r [^sgW + Ay5>M^''] 



(17) 



where F^u = di^uJv - d^cjfj, and Gfj.^ = dfj,p^ - d^p^ are 
strength tensors for uj field and p field, respectively, tp^ a, 
Ldfj,, p^ are nucleon field, isoscalar-scalar field, isoscalar- 
vector field and isovector-vector field, respectively, and 
the arrows denote the vector in isospin space. The Ag 
and Ay represent coupling constants between the isovec- 
tor p meson and the isoscalar a and lo mesons, respec- 
tively, which arc important for the description of the den- 
sity dependence of the symmetry energy. In addition, 
M is the nucleon mass and m^, m^j, mp are masses of 
mesons. 

In the mean field approximation, after neglecting ef- 
fects of fluctuation and correlation, meson flclds arc re- 
placed by their expectation values, i.e., a ^ ujq ~^ ^p,, 

-(3) -t 

Pq p^, where subscript "0" indicates zeroth compo- 
nent of the four-vector, superscript "(3)" indicates third 
component of the isospin. Furthermore, we also use in 
this work the non-sea approximation which neglects the 
effect due to negative energy states in the Dirac sea. The 
mean field equations are then expressed as 



m^cr =ga 



2 -(3) 



2 3/'f3)\^ — 

PS ~ baM {gao) - Co- (gcrff) + {opp}, 1 As5<tCT 

(18) 

P- (3a;Wo)^ - Ay5„(:jo (spPo^^) (19) 



Pp- Pn- As5pPo^^ ig<T(Ty - Av9pPo'' i9i^^o) 

(20) 



(3) 



K 



asy 



AT. 



sym 



6A, 



(16) 



and this expression has been extensively used in the lit- 
erature to characterize the isospin dependence of the in- 
compressibility of asymmetric nuclear matter in the lit- 
erature [1,11, HHIS^. 



where 

P = (V'7°V') = Pn + Pp, PS = (-ipip) = PSn + PSp, (21) 
are the baryon density and scalar density, respectively. 
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with the latter given by 



PSJ 



(2^)^ .0 



dk 



MX 



M*j 
2^ 



fc|2 + A/* 2 





J* 



, J = p, i(22) 



In the above expression, we have Ep 
with Mq — M ~ g^a being the nucleon Dirac mass and 
the Fermi momentum kp ~ kpll+r^S^^^ with ra = 1 for 
neutrons and T3 = —1 for protons, and kp = (3tt'^ 
being the Fermi momentum for symmetric nuclear mat- 
ter. 

The energy-momentum density tensor for the interact- 
ing Lagrangian density (|17p can be written as 

= 4,i-f^'^'''^lJ + ^^'a^''a-F^''^^''u;,J-G'''^^''p,J-Cg^"', 

(23) 

where g^^ = (+, — , — , — ) is the Minkovski metric. In the 
mean field approximation, the mean value of time (zero) 
component of the energy-momentum density tensor is the 
energy density of the nuclear matter system, i.e., 



2 -2 I 2 -2 I 2 1 -(3; 



113 



+ I {gpPo^y l^sigaa)^ + SAvig^Coof] , (24) 



where 



(27r) - JO 



dk^/\k\^ + Mf 



— I k^dk 



= -^{ZEi*pj^M*psj\, J = p,n, (25) 

is the kinetic part of the energy density. Similarly, 
the mean value of space components of the energy- 
momentum density tensor corresponds to the pressure 
of the system, i.e., 

1 

" 2 

- '^b^ig^af - ^c„{g„af + ^c^{g^u)of 

+ \ [gppf^ [As(3.^)' + Av(5c.cSo)'] , (26) 



2 -2 

ra„G 



,2-2 



,(3) 



where the kinetic part of pressure is given by 



1 
3^2 



dk 



k" 



fc2 + Mf 



I, J = p,n. 



(27) 



The binding energy per nucleon of the asymmetric nu- 
clear matter can be calculated through the energy density 

by 



Eip^5)^'-^-M. 



(28) 



Furthermore, the symmetry energy i?sym(p) can be ob- 
tained as 



Esymip) 



1 d^E{p,5) 



2! 9(52 



<5=0 



GEt 



g^pp 

2Qp 



(29) 



while the 4th-order symmetry energy Esym.iip) can be 
expressed as 



Esym,4{p) 



where 



1 d'^E{p, S) 



4! 



- Esym,4{p) + E^^m^Ap) 



(5=0 



648 



rp*5 



(30) 
(31) 



2g4 V Qu Q. 



glpM^kl ( 4AsgagpP _ M^kj, 
Ql 3Ef 



2AQaEf 



(32) 



with iJ^yJ^ 4(p) representing the kinetic part (including 
the interactions due to the nucleon effective mass) while 
■^sy-miip) the other part due to the interaction in the 
4th-order symmetry energy, and E*p = \JM^^ -f fc^. The 
coefficients Qcr, and Qp arc defined as 



Qo 



gl 



3ps 3p\ 3_ 
^ - ^ ) + 26.A/g.<T 



2,0 4-2 



=to2 _^ Ksglgla"^ + kyglgl^l- 



o 4-2 

(33) 
(34) 
(35) 



In the above expressions, all the fields are calculated in 
the case of symmetric nuclear matter, i.e., at (5 = 0. 

The analytical expression of the symmetry energy 
Esym{p), i-e., Eq. (f^ . is a well-known result firstly 
given in Ref. [38[. To our best knowledge, the formulas 
give, for the first time, the analytical expres- 
sion of the 4th-order symmetry energy Esym,4:{p) in the 
RMF model, which are the main results of the present 
work. These analytical expressions allow us to evalu- 
ate accurately the 4th-order symmetry energy Esym.iip) 
and thus study the higher-order corrections to the em- 
pirical parabolic approximation within the framework of 
the RMF model. Before presenting numerical results, it 
is instructive to analyze firstly the low density behavior 
of the EgymAip)- When p — >■ 0, the magnitude of all 
fields will approach to zero and both Ep and A/q will 



5 



approach to M, leading to E^^-^ ^{p) from Eq. ((32 



and £'sym.4(p) ^ AT ^^'om Eq. (PT|) . Therefore, m the 
low density limit, we have 



lim E. 



'sym,4 



1 



i-2 



162 M ^ 



(36) 



which is exactly the result from the free Fermi gas model 
as expected. 



III. RESULTS AND DISCUSSIONS 

A. The 4th-order symmetry energy and 
higher-order effects on the isobaric incompressibihty 
of asymmetric nuclear matter 
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FIG. 1: (Color online) Density dependence of the 4th- 
order symmetry energy -E'sym,4(p) as well as its kinetic part 
£'sym,4(p) and interacting part i5^m,4(p) from two accurately 
calibrated interactions, i.e., FSUGold (a) and lU-FSU (b). 

Shown in Fig. [T]is the density dependence of the 4th- 
order symmetry energy i?sym,4(p) as well as its kinetic 
part E^^jjyp) and interacting part E^^-^^^{p) using two 
accurately calibrated interactions, i.e., FSUGold [3y| and 
lU-FSU 0. The FSUGold has been accurately cali- 
brated to the ground-state properties of closed-shell nu- 
clei, their linear response, and the structure of neutron 
stars while the lU-FSU is a recently developed effective 
interaction that improves the FSUGold by incorporating 
some of the recent constraints on properties of neutron 
stars. One can see from Fig. [T]that the i?syin,4(p) is quite 
small (less than 0.7 MeV) at normal nuclear matter den- 
sity while it increases with density and can reach to about 
7 MeV at p = 1 fm~'^. Furthermore, one can see that 
the kinetic part £'^"^5^ ^{p) dominates over the interacting 
part iJ^jjj 4(/9) with the latter is generally negative. 

In order to investigate higher-order i?sym,4(p) effects on 
the EOS of asymmetric nuclear matter, it may make more 



FIG. 2: (Color online) Ratio of the 4th-order symmetry en- 
ergy to the symmetry energy as a function of density from two 
accurately calibrated interactions, i.e., FSUGold and lU-FSU. 



sense to calculate the ratio of the 4th-order symmetry en- 
ergy to the symmetry energy, i.e., Esy^^i{p) / Esy,^{p) . In 
Fig. [21 we show this ratio as a function of density with 
interactions FSUGold and lU-FSU. It is seen that the 
ratio i?sym,4 (p) /Esym [p] has a very small value of about 
2% around normal nuclear matter density po, but it can 
reach to about 6% '-^ 7% at high densities (e.g., 1.0 fm^'^ 
or 6 ~ 7po). This result is essentially consistent with 
the nomelativistic calculations in some phenomenologi- 
cal models [4l|. These features imply that the i?sym,4(p) 
may become important at higher densities in some ex- 
treme physical conditions such as in neutron star where 
the isospin asymmetry S can be close to unity. As an 
example, in the next subsection, we shall study effects 
of the 4th-order symmetry energy on the proton fraction 
in /3-stablc neutron star matter. In addition, one can 
see from Fig. [2] that in the low density limit, we have 
lim i?sym,4(p)/£'svm(p) = 1/27 as expected from the free 

Fermi gas model. 

The analytical expression of the 4th-order symmetry 
energy Egym.iip) allows us to calculate accurately the 
density slope and curvature parameters of i?sym,4(p), 
i.e., Lsyi„,4 and i^sym,4, and thus obtain the accu- 
rate value of the higher-order isobaric incompressibil- 
ity KsntA of asymmetric nuclear matter according to 
Eq. (Uni). Table U displays the characteristic parame- 
ters of asymmetric nuclear matter, namely, po, Eo{po), 

-E'sym(Po)j -E'sym,4(Po)i Jq: Iq ■ L, Kgym, Jsym: Isym: 

£sym,4j -K^sym,4i '^sym,4j Kasyi -K^sat,2j K^at^i and the ra- 
tios Ksa.t,2/Ki,sy and Ksi,t,4:/ Ksa.t,2, for the two accurately 
cahbrated interactions FSUGold and lU-FSU. To see the 
variation of the higher-order characteristic parameters 
with the density dependence of the symmetry energy, 
we also include in Table |T] the results from 5 interac- 
tions denoted as FSU-I, FSU-II, FSU-III, FSU-IV and 
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TABLE I: Characteristic parameters of asymmetric nuclear matter, namely, po (fm Eo{po) (MeV), i5sym(po) (MeV), 
£;sym,4(po) (MeV), Ko (MeV), Jo (MeV), lo (MeV), L (MeV), K,ym (MeV), Jsym (MeV), Jsym (MeV), L.ym,4 (MeV), ^-=^^,4 
(MeV), Jsym,4 (MeV), K^sy (MeV), iCsat,2 (MeV), /<'aat,4 (MeV), and the ratios Ks^t,2/K^y and K^^iaI Ks^t,2, for diflferent 
interactions. 





FSUGold 


lU-FSU 


FSU-I 


FSU-II 


FSU-III 


FSU-IV 


FSU-V 




0.148 


0.155 


0.148 


0.148 


0.148 


0.148 


0.148 


Eo{po) 


-16.3 


-16.4 


-16.3 


-16.3 


-16.3 


-16.3 


-16.3 


Esym{po) 


32.5 


31.3 


37.4 


35.5 


33.9 


31.4 


30.9 


-E'sym,4(po) 


0.66 


0.67 


0.66 


0.66 


0.66 


0.66 


0.78 


A-o 


229.2 


232.3 


229.2 


229.2 


229.2 


229.2 


229.2 


Jo 


-521.6 


-288.5 


-521.6 


-521.6 


-521.6 


-521.6 


-521.6 


Jo 


2815.7 


4541.9 


2815.7 


2815.7 


2815.7 


2815.7 


2815.7 


L 


60.4 


47.3 


109.5 


87.4 


71.7 


52.1 


49.4 


Ksym 


-51.4 


29.0 


2.7 


-68.4 


-74.4 


-16.7 


5.5 


J 3ym 


426.5 


363.9 


-101.4 


157.6 


399.0 


251.2 


80.4 


Jsym 


6331.8 


11346.5 


285.8 


1364.1 


4211.9 


6136.2 


4620.4 


J/sym,4 


1.9 


1.8 


1.9 


1.9 


1.9 


1.9 


2.3 


J^syni,4 


0.5 


0.1 


0.5 


0.5 


0.5 


0.5 


0.1 


Jsym ,4 


5.0 


6.3 


4.8 


4.8 


4.8 


5.2 


5.1 


J ^ asy 


-413.8 


-255.1 


-654.3 


-592.6 


-504.4 


-329.4 


-290.9 


Ji'sat,2 


-276.4 


-196.3 


-405.1 


-393.8 


-341.4 


-210.8 


-178.5 


Ji'sat,4 


3.0 


47.6 


338.6 


183.4 


50.0 


12.9 


32.4 


Ji^sat,2/A'asy 


67% 


77% 


62% 


66% 


68% 


64% 


61% 


— KbuX^a/ JS^sat,2 


1% 


24% 


84% 


47% 


15% 


6% 


18% 



FSU-V for which the parameters {As, Ay) are selected 
as (0.00,0.00), (0.00,0.01), (0.00,0.02), (0.00,0.04) and 
(0.01,0.03), respectively, while the gp parameter is ad- 
justed accordingly to fix EgymiPf) = 25.57MeV at pf = 
0.1 fm""^ as in the FSUGold interaction. The other pa- 
rameters for FSU-I, FSU-II, FSU-III, FSU-IV and FSU-V 
are exactly the same as in FSUGold (Note: the FSUG- 
old corresponds to the case of (As, Ay) = (0.00,0.03)). 
This is equivalent to solve a constraint equation about 

2 

gp, As and Ay, i.e., = 13.47MeV where a and loq in 
Qp (See Eq. f ([55|) )) are determined by the properties of 
symmetric nuclear matter in FSUGold. 

From Table m one can see that the i?sym,4(Po) is gen- 
erally less than 1 MeV (about 0.66 MeV for most of the 
interactions considered here), consistent with that ob- 
served in Fig. [1] These results about the i?sym,4(po) are 
further in agreement with the calculations in the nonrel- 
ativistic models of MDI and Skyrme-Hartree-Fock j28| . 
nicely verifying the empirical parabolic law around the 
normal nuclear matter density po- 

As pointed out previously, the difference between 
Ksat,2 and J^asy reflects the contribution from higher- 
order effects, namely, the value of JqL/Kq, which has 
been usually neglected in many calculations in the lit- 
erature 0, i, [il [H, Hi]. From Table HI one can see 
that neglecting the JqL/Kq term generally leads to 30- 
40% relative error for the JCsat,2 parameter and thus the 
JqL/Kq term contribution to the ii'sat,2 parameter can- 



not be neglected simply, confirming the previous findings 
in the nonrelativistic studies [28j . 

Furthermore, it is seen from Table U that the value of 
higher-order Jirsat,4 is generally small compared with that 
of Ksat,2 for most of the interactions considered here. 
In addition, one can see that the i^sat,4 becomes more 
important for the interactions with larger L values and 
this is consistent with the nonrelativistic studies [2^. It 
should be noted that the higher-order ifsat,4 term can 
be safely neglected in the study of giant resonance of 
finite nuclei [42| where the isospin asymmetry S is usually 
small, i.e., about 0.2. 

B. Effects of i5sym.4 on the proton fraction in 
/^-stable nuclear matter 

In order to further illustrate the effects of the 4th- 
order symmetry energy on the EOS of asymmetric nu- 
clear matter, we calculate the proton fraction Xp in /3- 
stable neutron star matter where the isospin asymmetry 
6 is generally close to 1. The chemical composition of the 
neutron star is determined by the requirement of charge 
neutrality and equilibrium with respect to the weak in- 
teraction (/3-stable matter). From the binding energy per 
nucleon, i.e., Eq. ([1]), we can calculate the proton frac- 
tion, Xp = (1 — (5)/2, for /3-stable nuclear matter as found 
in interior of neutron stars. For neutrino free /3-stable nu- 
clear matter, the chemical equilibrium for the reactions 
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n — !■ p + e + Ve and p 



e n + Ve requires 

dE 
'dS 



= UE,yn,{p) + SS^'E.yn.^^ip) + 0{S^) (37) 

where /x^ = dEi/dxi [i = n,p,e,fi) is the chemical po- 
tentiaL For relativistic degenerate electrons, we have 



k 



ml + (37rVa;e)'/^ 



1/2 



,1/3 



(38) 



where me = 0.511 MeV is the electron mass, and Xp = Xe 
because of charge neutrality. 

Just above a nuclear matter density at which ex- 



= 0.105 GeV, the reactions 
n + u^ and n 



ceeds the muon mass m 

e~— >-/i~+z/e + j^^,p + /x" n + v^ ana n + 
are energetically allowed so that both electrons and 
muons are present in /3-stablc nuclear matter, this alters 
/3-stability condition to 




FIG. 3: (Color online) Density dependence of the proton frac- 
tion Xp in /3-stable npe/i matter with FSUGold (a) and lU- 
FSU (b). Three cases, i.e., the full EOS of asymmetric nuclear 
matter (solid lines), its parabolic approximation (up to 5^ in 
Eq. ID)) (dotted lines), and further including the 4th-order 
symmetry energy (up to 5* in Eq. ([1])) (dashed lines), are 
considered. 

In Fig. [3l we show the density dependence of the pro- 
ton fraction Xp in /3-stable npep, matter with the inter- 
actions FSUGold and lU-FSU. We consider three cases 
for the EOS of asymmetric nuclear matter here, i.e., the 
full one. its parabolic approximation (up to 5"^ in Eq. 
([1])), and the one further including the 4th-ordcr sym- 
metry energy (up to 5'^ in Eq. HI)). The results show 
that, for both interactions of FSUGold and lU-FSU, the 



4th-ordcr symmetry energy is moderately important for 
the proton fraction, especially at higher densities. For 
the FSUGold (lU-FSU) interaction and p ^ 1.0fm~^ 
for instance, including the 4th-order symmetry energy in 
the parabolic approximation to the EOS of asymmetric 
nuclear matter will increase the proton fraction Xp from 
15.37% (13.44%) to 16.43% (14.81%), producing a rela- 
tive variation of about 7% (10%). These results indicate 
that the 4th-order symmetry energy may have obvious 
effects on the proton fraction Xp in /3-stablc npep, matter 
and the parabolic approximation to the EOS of asym- 
metric nuclear matter may significantly underestimate 
the proton fraction, especially at higher densities. These 
features are consistent with the nonrelativistic Skyrme- 
Hartree-Fock calculations [25l |. 

Furthermore, one can see from Fig. [3] that the dif- 
ference between the results with the full EOS and with 
the one containing the terms up to the 4th-ordcr sym- 
metry energy is very small, indicating that the EOS of 
asymmetric nuclear matter including the terms up to the 
4th-order symmetry energy (up to 5^ in Eq. ([1])) could 
be a good approximation for the determination of the 
proton fraction in /3-stable npefi matter. 



C. Effects of i?sym,4 On core-crust transition density 
and pressure in neutron stars 

The transition density pt is the baryon number density 
that separates the liquid core from the inner crust in neu- 
tron stars and it plays an important role in determining 
many properties of neutron stars [13, [H, |43| - [46| . One 
simple and widely used way to determine the core-crust 
transition density pt is the so-called thermodynamical 
method, which requires the system to obey the following 
intrinsic stability condition |46l - l48| 



dP 
dv 

df-L 



np 



> 0, 



> 0, 



(40) 
(41) 



where the P = Pi, + is the total pressure of the npe 
matter system with Pb and P^ denoting the contributions 
from baryons and electrons respectively, and the v and 
Qc are the volume and charge per baryon number. The 
finp is defined as the chemical potential difference be- 
tween neutrons and protons, i.e., pnp = Pn — Hp- The 
pressure Pe is only a function of the chemical poten- 
tial difference pnp by assuming the /3-equilibrium con- 
dition is satisfied, i.e., pnp = Me- By using the relation 
dEb{p,Xp)/dxp = -finp with Eb{p,Xp) being energy per 
baryon from the baryons in the /3-equilibrium neutron 
star matter and Xp ~ Pp/ p, and treating the electrons 
as free Fermi gas, one can show ^3] that the thermo- 
dynamical relations Eq. (PU)) and Eq. (|1T]) are actually 
equivalent to the following condition 
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thermal 



2p 



dEb{p,Xp) 2d'^Eb{p,Xp) fd'^Ei,{p,Xp 



dp 



dp^ 



dpdxp 




d'^E(,{p,Xp) 
dx'^ 



>0, 



(42) 



which determines the thermodynamical instabiUty region of the /3-equihbriuni neutron star matter. The baryon 
number density that violates the condition Eq. (j42p then corresponds to the core-crust transition density in neutron 
stars for the thermodynamical method. 

With the EOS of asymmetric nuclear matter including the terms up to the 4th-ordcr symmetry energy (i.e., up to 
(5"* in Eq. ([T])), Eq. (g^]) is then reduced to 



Vt 



thermal 



, d^Eo 
dp^ 



2p- 



dEn 



dp 



p^S^ 



2 d'^Esyn-i{p) 

' dp^ 

dEsymip) 



dEsymip) 



Esym{p) + 6i?sym,4(p) 



dp 



26 



Op 

2 dEsyin,4{p) 



dp 



5^ 



> 0. 



> d"^ EsymA{p) , r,^9i?sym,4(p) 



dp^ 



+ 2/9- 



dp 



(43) 



The baryon number density that violates the condition Eq. (j43p then corresponds to the core-crust transition density 
in neutron stars for the EOS of asymmetric nuclear matter including the terms up to the 4th-order symmetry 
energy (up to 5^ in Eq. ([1])). The corresponding transition pressure P/*^ at p^^^ for the EOS of asymmetric nuclear 
matter including the terms up to the 4th-order symmetry energy (up to 5^ in Eq. ([l}) is then given by 



-,fc,4th 



1 ^e,4th 










^-^sym 




dEsyn\,4 


dp 


dp 


dp 


OEo^^, 


dEgym 


f<54 


dEgyiT^4 


op 


dp 


dp 



p=pt 



p=pt 



ApS 



1-S 



• [Esym + '^^£'sym,4] 



(44) 



p=pt,S=St 



where St is the isospin asymmetry of the /3-cquilibrium neutron star matter at the corresponding transition density. 

The transition density obtained by neglecting the -Esym. 4 (p) term in the condition Eq. (|43p corresponds to the 
core-crust transition density Pj"'^ in neutron stars for the parabolic approximation to the EOS of asymmetric nuclear 
matter (up to 6^ in Eq. ([T])). The corresponding transition pressure at p^nd -j-qj. ^-j^^ parabolic approximation to the 
EOS of asymmetric nuclear matter (up to in Eq. ([1])) is then expressed as 



p2nd _ 



pb,2nd . pe,2nd 



dEo 
dp 



j2 dEsyra 

dp 



pf'^Pe 



P=Pt 



[ dEo 
\ dp 



^2 OEgyi^ 

dp 



p=pt 



ApS 



2 ■ J^sym 



(45) 



p=Pf,5=St 



Due to simplicity, the p^'^'^ and P^nd j-^g^yj, been exten- 
sively applied to determine the inner edge of neutron star 
crusts within the nonrelativistic models [i^, l48l - [50| and 
recently in the RMF model [5l| as well. However, re- 
cent studies based on some nonrelativistic models have 
demonstrated [27| that the parabolic approximation to 
the EOS of asymmetric nuclear matter may lead sys- 
tematically to significantly higher core-crust transition 
densities and pressures, especially with stiffer symmetry 
energy functionals. It is thus very interesting to see how 
the higher-order Esym.iip) affects the transition density 
Pt and pressure Pt in the RMF model. 

In Table HI we show the p?"^, pf^, P^""^ , and ^ ob- 
tained from the thermodynamical method with different 
interactions as in In Table HI It is interesting to see that 
including the 4th-ordcr symmetry energy in the parabolic 



approximation to the EOS of asymmetric nuclear matter 
indeed reduces significantly the core-crust transition den- 
sity Pt , which is consistent with the nonrelativistic calcu- 
lations [l^]- Furthermore, one can see that the 4th-order 
symmetry energy may have even more drastic effects on 
the core-crust transition pressure Pt, namely, including 
the 4th-order symmetry energy in the parabolic approxi- 
mation to the EOS of asymmetric nuclear matter reduces 
drastically the core-crust transition pressure Pt- There- 
fore, our results indicate that the empirical parabolic ap- 
proximation may cause large errors for the determination 
of the Pt and Pt in neutron stars in the nonlinear RMF 
model. 



9 



TABLE II: The pf"^ (fm-=*), ' (fm"^), Pt""^ (MeV/fm^), and ' (MeV/fm''') obtained from the thermodynamical method 
with different interactions. 





FSUGold 


lU-FSU 


FSU-I 


FSU-II 


FSU-III 


FSU-IV 


FSU-V 




0.089 


0.090 


0.085 


0.088 


0.088 


0.083 


0.080 




0.051 


0.077 


0.069 


0.054 


0.053 


0.068 


0.072 


p2nd 


1.316 


0.673 


0.664 


1.010 


0.968 


0.621 


0.501 


p4th 


0.321 


0.530 


0.302 


0.236 


0.259 


0.420 


0.414 



IV. SUMMARY 

We have derived for the first time the analytical expres- 
sion of the nuclear matter fourth-order symmetry energy 
Esym,A{p) within the framework of the nonlinear RMF 
model. It should be mentioned that the analytical ex- 
pression of Esyin,iip) can be easily generalized to the case 
of the density dependent RMF model that has similar 
isospin structure as the nonlinear RMF model (See, e.g., 
Ref . [Ill ) . This provides the possibility to investigate the 
higher-order -Esym,4(p) corrections to the widely used em- 
pirical parabolic law for the asymmetric nuclear matter 
in the RMF model. In the present work, as examples, 
we have investigated the Esym,4{p) effects on the prop- 
erties of asymmetric nuclear matter, the proton fraction 
Xp in /3-stablc npefi matter and the core-crust transition 
density pt and pressure Pt in neutron stars within the 
nonlinear RMF model with two accurately calibrated in- 
teractions, i.e., FSUGold and lU-FSU. 

Firstly, our results have indicated that the value of 
Esym,iip) at normal nuclear matter density po is gener- 
ally less than 1 McV, and thus the empirical parabolic 
approximation E{p,5) ~ Eo{p) + i?sym(p)<5^ has been 
nicely confirmed around pQ. However, at higher densi- 
ties such as 1 fm~^, the value of Eeym^ip) can be about 
7 MeV and the ratio of i?sym,4(p)/£'sym(p) can reach to 
about 7%. These results imply that the i?sym,4(p) may 
become nonnegligible at higher densities. Furthermore, 
the analytical form of the i?sym,4(p) allows us to study 
the higher-order effects on the isobaric incompressibility 
of asymmetric nuclear matter. Our results have indi- 
cated that the value of higher-order Ksat,4 is generally 
small compared with that of ifsat^ confirming the pre- 
vious nonrelativistic calculations |28|. 

Secondly, for the proton fraction Xp in /3-stable npep, 
matter, we have found that, compared with the results 
from the empirical parabolic approximation to the EOS 
of asymmetric nuclear matter, including the 4th-ordcr 
symmetry energy i?sym,4(p) can enhance the proton frac- 
tion Xp by about 10% at higher densities. These results 



indicate that the empirical parabolic approximation to 
the EOS of asymmetric nuclear matter may cause obvi- 
ous errors for the determination of the proton fraction 
in neutron stars within the nonlinear RMF model, which 
is in agreement with the results from the nonrelativistic 
models [25|. 

Finally, we have demonstrated that including the 4th- 
order symmetry energy £'sym.4(/o) in the parabolic ap- 
proximation to the EOS of asymmetric nuclear matter 
can reduce significantly the core-crust transition density 
Pt and furthermore it has even more drastic effects on 
the core-crust transition pressure Pt- Therefore, our re- 
sults have clearly demonstrated that the extensively used 
empirical parabolic approximation to the EOS of asym- 
metric nuclear matter may lead systematically to signif- 
icantly higher core-crust transition density pt and pres- 
sure Pt in neutron stars within the nonlinear relativistic 
mean field model, confirming the previous finding based 
on nonrelativistic calculations 



27 



Therefore, we conclude that the higher-order Esym.iip) 
in the EOS of asymmetric nuclear matter may have dif- 
ferent effects on different quantities, and generally one 
cannot simply neglect them, especially under some ex- 
treme physical conditions, such as in neutron stars. 
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